Emergent spin excitations in a Bethe lattice at percolation 
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We study the spin 1/2 quantum Heisenberg antiferromagnet on a Bethe lattice diluted to the percolation threshold. Dilution creates areas of 
even/odd sublattice imbalance resulting in "dangling spins" (L. Wang and A. W. Sandvik, Phys. Rev. Lett. 97, 117204 (2006); Phys. Rev. 
B 81, 054417 (2010)). These collectively act as "emergent" spin 1/2 degrees of freedom and are responsible for the creation of a set of low 
lying "quasi degenerate states". Using Density Matrix Renormalization Group (DMRG) calculations, we detect the presence and location of 
these emergent spins. We find an effective Hamiltonian of these emergent spins, with Heisenberg interactions that decay exponentially with 
the distance between them. 



Introduction — Quantum spins on percolation clusters [1] 
provide an ideal testbed for studying the interplay between 
geometrical disorder and quantum fluctuations. The Hamilto- 
nian for these problems is 

where Si are Pauli spin 1/2 operators and the sum runs over 
nearest-neighbor occupied sites, and J > 0. Theoretical [2, 3] 
and experimental [4] studies of quantum spins on diluted 
square lattices have focussed on the question of whether long 
range order survives up to the classical percolation threshold 
Pc- A numerical study [5] seems to have settled this question 
and found long range order to be robust to quantum fluctua- 
tions, surviving all the way up to pc- 

The excitations are less straightforward. For uniform lat- 
tices with number of sites N, the lowest energy scale con- 
sistent with Neel order breaking a continuous symmetry is 
~ JN~^, corresponding to a "tower" of states: mixtures 
of symmetry-broken states that become degenerate in the 
thermodynamic limit [6-8]. However a Quantum Monte 
Carlo study by Wang and Sandvik [9] discovered a some- 
what "anomalous" finite size scaling of the spin gap A/o^,: 
A/oui ~ A^^^ (for clusters with a singlet ground state) or 
A/oM) ~ N^^-^ (for generic clusters, most with a non-singlet 
ground state). A strong coiTespondence was shown [9] be- 
tween these low lying states and places on the cluster where 
there is a local imbalance between the number of even and 
odd sites. It was conjectured that, in each such place, a spin 
degree emerges which is effectively decoupled from the anti- 
ferromagnetic order and hence was called a "dangling spin." 

The goal of this Letter is to characterize the dangling-spin 
degrees of freedom numerically, relating their nature to the 
local geometry of the cluster, and to explain the observed low 
energy spectrum in terms of mediated interactions between 
dangling spins. Our Hamiltonian is (1) on clusters obtained 
by randomly diluting the Bethe Lattice of coordination 3 at its 
percolation threshold, pc — 1/2 (see examples of small clus- 
ters in Fig. 1). The lack of loops in the Bethe lattice is con- 
ducive for using the Density Matrix Renormalization Group 
(DMRG) algorithm [10], as adapted to generic tree graphs, to 
obtain ground and (some) excited states [11]. 

In the rest of this Letter, we first show that a typical per- 
colation cluster's spectrum has a clearly separated low energy 
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Figure 1: (Color online) Three different percolation clusters (all of 
the same size — 18) are shown with their corresponding low en- 
ergy spectra. The 'red' and 'green' circles indicate even and odd 
sites. The broken dashed lines show dimer coverings which serve 
as a heuristic to locate the "dangling spins". Dangling spins in each 
cluster have been circled with thick black lines. Low energy spec- 
tra for each of the clusters show low lying quasi degenerate (QD) 
states separated from the continuum by an energy scale A. ago (not 
shown) is a measure of the spread of energies in the QD space. 

component, with a multiplicity consistent with the expected 
number of weakly coupled spin-1/2 (sometimes spin-1) dan- 
gling spins. We next show that each dangling spin is some- 
what delocalized over a few sites: on the one hand, we model 
it as an unpaired spin in a dimerized background to predict the 
dangling spin's nature from the local geometry; on the other 
hand, by processing spin expectations we obtain the explicit 
"localized state wavefunction" for each dangling spin. Finally, 
for each cluster we construct the effective Hamiltonian of the 
emergent dangling spins, consisting of pairwise, unfrustrated 
exchange interactions decaying exponentially with separation, 
mediated by the background of almost dimerized spins on the 
balanced parts of the cluster; this accurately reproduces the 
details of that cluster's low-energy spectrum. 

Exact correspondence between dangling spins and low en- 
ergy spectrum — We caiTied out DMRG calculations for sev- 
eral hundred balanced clusters (i.e. having equal number of 
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even and odd sublattice sites) [12] for sizes up to iV = 100, 
targeting multiple excited states in the low energy spectrum. 
Since the number of low energy states was found to increase 
rapidly with an increase in the number of dangling spins, we 
restricted our analysis to the case of four dangling spins [13]. 
For analyses where only the ground state properties were of 
interest or where only the lowest energy gap Ajo^ was re- 
quired, larger clusters {N = 200), with an arbitrary number 
of dangling spins, were also studied. 

In a typical percolation cluster, we observed a distinct set of 
low-lying energy levels we shall call "quasi-degenerate" (QD) 
[14] since (we claim) they would be exactly degenerate in 
the limit that the dangling spins are fully decoupled from the 
rest of the sites. The QD states are separated from the contin- 
uum of higher energy Quantum Rotor and spin wave states by 
a finite size gap we call A (specifically defined as the differ- 
ence between the mean of QD levels and the lowest non-QD 
level). The set of QD states are identified by looking at the 
difference in energies of consecutive states up to the Quantum 
rotor excitation and finding the pair of states with the largest 
gap. The lower energy state in this pair and all states below 
that make up the QD spectrum [15]. The energy scale charac- 
terizing the spread of the QD states, ctqd, is defined to be the 
standard deviation of the QD energies from their mean value. 
The ratio r = ctqd/A was found to be small (for example 
r ~ 0.17 ± 0.1 for N — 50), justifying our notion of a sepa- 
ration of scales. 

Figure 1 also shows a striking correspondence between 
the number of low lying QD states, A/qd. and the number 
of dangling spins Ud on the percolation cluster We find 
that A/qd = 2"i/2 3"\ where ni/2 and ?ii are integers and 
ni/2 + 2ni = Ud- Our interpretation of this multiplicity is 
that 2ni of the dangling spins pair up so as to form a spin- 
1, while the others remain as spin-half degrees of freedom. 
There is thus a one-to-one correspondence between the low- 
energy (QD) eigenstates and the Hilbert space of the posited 
emergent spins. We used an algorithm (to be described in the 
next section) that relies only on the cluster geometry to objec- 
tively predict the numbers ni/2 and ni for each cluster, and 
verified that their predicted relationship with A/qd was satis- 
fied in every cluster. 

We also directly measured the lowest singlet-triplet gaps 
Alow for an ensemble of balanced clusters. Its typical value 
scales as ]\[-^-^±o.i^ which appears remarkably similar to the 
scaling previously seen on square lattice percolation clusters 
[9]. 

Locating Dangling degrees of freedom in real space — Hav- 
ing established the presence of emergent spin-half and spin- 
one degrees of freedom, we now develop two complementary 
ways of looking at them. 

The first is within the framework of a quantum monomer- 
dimer model. We imagine that the wavefunction is a product 
of valence bonds in which the N spins are paired (dimerized) 
into singlets to the maximum extent possible (optimal config- 
uration). Even when even and odd sites are balanced globally, 
there remain some uncovered sites, i.e. monomers, due to lo- 
cal imbalances. These are spin- 1/2 degrees of freedom and 
(within this picture) represent the dangling spins. There are 



multiple ways to optimally place the monomers; the actual 
wavefunction is imagined to be a superposition of these ways. 

Our geometric algorithm, based on the valence bond frame- 
work, finds one element from the set of optimal dimerizations 
of the cluster and then attempts to find other elements of the 
set by locally replacing monomers with adjacent dimers. In 
spirit, this is a "greedy" algorithm which tries to place dimers 
wherever possible (to obtain an optimal dimerization pattern), 
working from the outer sites inwards on the cluster. 




Figure 2: (Color online) :Typical geometrical motifs in Cay ley tree 
percolation clusters, as related to monomer/dimer construction. We 
disconnect the cluster into "spans", shown by (blue) dashed circles, 
at the "prong" bonds, as indicated by the (red) cut lines. The thick- 
ness of the (grey) bonds is directly proportional to {Si ■ Sj). 




Figure 3: (Color online): Spatial profiles associated with "dangling 
spins", shown on the subclusters marked (a)-(d) of the same cluster 
shown in Fig. 2. The area of the black circles is proportional to ua 
given by Eq. (4). (a) shows a "fork" (b) shows a site surrounded by 
two "hairs" (c) shows a subcluster where two monomers on the same 
sublattice are present forming an effective spin- 1 . The (blue) arrows 
indicate the monomer is free to hop around (delocalize) within the 
subcluster. (d) shows a region where the spins are "inert" (largely 
dimerized) and do not participate in the low energy physics. 

Given any cluster, there are two operations which cut it 
down to a smaller cluster or clusters, such that all optimal 
dimerizations on the smaller cluster(s) are in 1-to-l correspon- 
dence with some of dimerizations on the larger one. The first 
operation is that, wherever two sites have coordination 1 or 
2, we can remove both (given the dimerization on the smaller 
cluster, just insert another dimer to get the dimerization on 
the larger one). The second operation is that, wherever we 
find a pair of adjacent sites with respective coordinations 3 
and 1 (a "prong"), we can always place a dimer on that pair. 
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which fragments the rest into two subclusters (Fig. 2); a very 
common special case is the fork, at which we can arbitrarily 
choose either side to be the "prong". These two operations 
can be used recursively till only isolated sites remain, each 
corresponding to one monomer in the original cluster Fur- 
thermore, any other optimal dimerization is accessible from 
the special one we constructed, by repeatedly exchanging a 
monomer with an adjacent dimer. 

A monomer can thus "hop" to sites on the lattice via such lo- 
cal monomer-dimer rearrangements as shown in Fig. 3(c). Our 
rule-of-thumb is that two monomers (of the same sublattice) 
form a spin-1 if and only if they can approach to the minimal 
separation of two steps [16]. We check this possibility by enu- 
merating all dimerizations accessible from the one generated 
by our algorithm. 

Fig. 2 also shows some typical geometrical motifs seen in 
an ensemble of percolation clusters. In Fig. 2(a), a spin-half 
dangling spin is localized at the tip of the fork. The arms of 
the fork have a strong dimerization and prevent the localized 
monomer from hopping to other sites in the cluster. Fig. 2(b) 
shows another type of spin-half degree of freedom, which is 
trapped in the interval between (or on) two "prongs". Fig. 2(c) 
is an example of an emergent spin-one excitation. For compar- 
ison, we also show in Fig. 2(d) a region of strong dimerization 
with no dangling spins. 

Our second way to capture the spatial nature of a dangling 
spin degree of freedom starts from the idea that it could be adi- 
abatically connected to an uncoupled spin, analogous to the 
Landau quasiparticle that is adiabatically connected to a free 
electron. Thus, our program is to label each emergent spin- 
1/2 degree of freedom by a "quasi spin" operator T^, where 
the index a labels each region on the cluster with a local spin 
imbalance. The T^'s are idealized as having a spin-half alge- 
bra. The actual "quasi spin" excitation is a composite object 
involving multi-spin excitations, localized on a few sites. 

Our assumption is that the quasi-spin quantum numbers are 
sufficient to label all the QD states; furthermore, we expect 
the action of any spin operator S^, when restricted to the QD 
states, practically reduces to a linear combination of Tq,'s act- 
ing on the quasi spins. Specifically, let Pqd be the projection 

operator onto the QD subspace. Let S^ qd = ^QDSi-PqD 
Then, 
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where each mode has most of its weight on sites within 
the region a and is expected to decay rapidly outside [17, 18]. 

Two operators P, Q are said to be orthogonal when their 
Frobenius inner product (^P,Qj = Tr ^P^q) is exactly 
zero. In this sense, the operators are orthogonal to each 
other. Since each is a quasi-spin 1/2 operator its inner 
product with itself is 1/2 (for each spin component). 

In light of Eq.(2), we can also construct a good approxima- 
tion Tq, to each operator T^, by choosing any representative 
site i in the dangling region a and normalizing the restriction 



of its spin operator to the QD states: 



(3) 



where ||0|| = (0,0)^ is the norm of any operator O. Note 
that the Tq,'s are not orthogonal to each other. A procedure to 
construct the Tq,'s from the T^'s will be discussed in the next 
section. 

Given the proposed relationship of the bare spins to the 

quasi spins, we discuss two related but independent measure- 
fa) 

ments to recover the mode vectors from numerically eval- 
uated expectations. For all the measurements proposed here, 
it was sufficient to record (within DMRG) all matrix elements 
of the form {l'\S^\l) and energies Ei for eigenstates 1,1' in 
the QD subspace. 

First, we consider the operator overlap ujij between two 
spins / and j on the lattice, defined to be. 
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(4) 



We substitute our ansatz (2) into (4) and use the operator or- 
thogonality of the Tq's, to get cuij — "j"''- 
consider a site / to be well within a dangling region a (i.e. 

(a) 

u\ is relatively large) then the amplitude on the remaining 
sites j (but far away from other dangling regions) is approx- 
imately Thus, the relative amplitudes of the mode 
vector can be recovered by this method. 

Our second measurement involves computation of the inter- 
site spin susceptibility matrix. 



Xij 



(Olefin) (n|^||0) 

En — Eq 



(5) 

where r is imaginary time, |0) denotes the ground state and 
En is the energy of an excited state \n) [19]. 

Though the sum runs over all excited states, it can be well 
approximated by taking only the states in the QD subspace. 
Then Xij can also be expressed in terms of the mode profiles 



Xij 



a/3 



X. 



(0|T^|n)(7i|TI|0) 



^eQD 



En — Ef) 



(6) 

Consider site in dangling region a (/3). From equation 
(6) it follows that Xij ~ uf^ u]p Xap, where the last factor 
is independent of sites i, j (so long as we stay within those 
regions). Within this approximation, the susceptibility matrix 
breaks up into blocks of rank 1 from which we can immedi- 
ately pull out the uf^ and uf'^ modes. 

Effective Hamiltonian in the Quasi degenerate subspace — 
According to our ansatz (2), there is a one-to-one correspon- 
dence between the QD Hilbert space and the Hilbert space 
of a set of abstract "quasi spin" operators T^. (For simplic- 
ity, assume they all have spin 1/2.) The latter are labeled us- 
ing an Ising basis \(j)t), where t stands for the quantum num- 
bers {if,t|, ...i^^}, with ta = ±1/2. We want to find the 
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unitary matrix M of coefficients expressing the QD states 
1^) (in eigenenergy basis) in terms of the quasi-spin basis, 

Using T„ from (3), we define Q^^'^ = (| ± T^), which 
is almost a projection operator, and let 

^ "d 

\^t)^X{Q'^\n) (7) 

Q = l 

where |n) could be any QD state (that is not annihilated by the 
operator prefactors) and \4>t) is normalized. Finally, define a 
matrix by Vltf = {4>t\4>t') - which is almost the identity 
matrix - and construct the orthonormal quasi spin basis of the 
QD states as 

\<t^t)^Y.^^~"^)u'\M. (8) 
t' 

The quasi spin operators Tq are then defined so as to have the 
appropriate matrix elements in this basis. 

Now consider the effective low energy Hamiltonian written 
in terms of the many body eigenstates 

;eQD «' 

where Ei is the eigenenergy of QD state |/), and the matrix el- 
ements htt' can be calculated since we know the transforma- 
tion between the bases {|^)} and {|(/)t)}. Every term \4>t){4't'\ 
can be uniquely expressed as a polynomial in the spin opera- 
tors {T„-}and{r±}. 

The effective Hamiltonian (9) then takes the new form of, 

Hes = ^ J^ivT^ti ■ T^v + multi spin terms; (10) 

(The two-spin terms must have this Heisenberg form due to 
the exact rotational symmetry retained by the QD states.) 

Although the magnitude of J^^ depends on the detailed ge- 
ometry of the cluster along the path connecting dangling re- 
gions fi and V, roughly speaking it decays with distance (using 
as metric the number of steps within the Bethe lattice, the so- 
called "chemical" distance). This is quantified by the scatter 
plots in Fig. 4, for an ensemble limited to clusters of equal 
size TV = 100 [20] each having two dangling spin- 1/2 spins 
{ud = ni/2 = 2). Since each dangling region is in general 
spread out over multiple sites, we must define an "effective 
distance" between two of them: 

where dij is the distance between sites i and j belonging 

to dangling regions and u, and the amplitudes uf^ for 
mode a are normalized. Figure 4 shows that indeed J^^ w 

Joe"'^"-/«, where (Jo,0 ~ (+0.15,4.85) for an even/odd 
pair of dangling spins, which are always antiferromagnetically 
coupled, or (—0.22, 4.55) for a pair on the same sublattice. In 
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Figure 4: (Color online) Effective couplings J12 between two dan- 
gling spins as a function of their effective separation obtained from 
clusters having exactly two dangling spins, (a) dangling spins on op- 
posite sublattices (antiferromagnetic coupling) (b) on same sublattice 
(ferromagnetic coupling). 

tht ferromagnetic case, choosing to fit only to clusters which 
do not form a spin-1 gives parameter values closer to the an- 
tiferromagnetic case. 

We considered another ensemble (not plotted) of clusters 
with = 50 having four dangling spins (n^ — 4) and 
obtained the effective Hamiltonian using the same prescrip- 
tion. We found that the non-pair terms typically account for 
a weight of at most 5% (using the Frobenius norm) in the 
Hamiltonian (10), confirming that the effective Hamiltonian is 
well approximated by pairwise Heisenberg exchange (at least 
in the limit of dilute monomer concentration). 

Conclusion — We addressed the excitations of the spin 1/2 
Heisenberg antiferromagnet on Bethe lattice percolation clus- 
ters. The emergent low energy degrees of freedom were found 
to be composite "dangling spin" excitations arising wherever 
there is a local imbalance between the even and odd sublat- 
tices [9]. (The same kind of imbalance was also recently 
shown to determine the low energy spectra of regular Cayley 
trees [11].) We classified dangling spins based on the local 
cluster geometry and described an algorithm to predict their 
locations on clusters of arbitrary size. 

Each of these emergent degrees of freedom is associated 
with a profile (m|"^ in the text) that plays the role of a "spinon 
wavepacket wavefunction" [21]; the degree of localization of 
the dangling spins could, in future work, be quantified using 
metrics (e.g. localization lengths) borrowed from Anderson 
localization. We did not, however, address the fundamental 
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reason why a dangling spin degree of freedom decouples from 
the cluster, an effectively zero-energy excitation; that will be 
addressed in a subsequent publication [17]. 

The dangling spins interact via an effective low energy 
Hamiltonian, consisting of unfrustrated Heisenberg couplings 
that decay (roughly) exponentially with separation. If we 
adopt the fitted dependences in Figure 4 as our definition of 
our Hamiltonian, it should be possible to understand the low- 
energy states for arbitrary clusters that are much too large for 
the DMRG (involving thousands of sites), but will be tractable 
using the strong disorder renormalization group method [22], 
and thus finally explain the scaling of the spin gap Aiow with 
cluster size that was found in Refs. 9. It would be interest- 
ing to see if the finite size scaling is exactly iV~^ and if so. 



whether it is universal. 

Finally, how does long range magnetic order survive in 
the scenario of dangling spins mediating interactions across 
strongly dimerized inert regions? We observed that spin corre- 
lations had an exponential decay across dimers but revived in 
regions on the cluster where a dangling spin was located. This 
suggests that in the thermodynamic limit, locally unbalanced 
regions are crucial for the propagation of magnetic order on 
percolation clusters. 
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